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, ^ . Abstract 

Masses of the pseudoscalar (J^*^ = ^) charmonium and bottomonium hybrids are determined using QCD 
' Laplace sum-rules. The effects of the dimension-six gluon condensate are included in our analysis and result in a 

, stable sum-rule analysis, whereas previous studies of these states were unable to optimize mass predictions. The 

pseudoscalar charmonium hybrid is predicted to have a mass of approximately 3.8 GeV and the corresponding 
bottomonium prediction is 10.6 GeV. Calculating the full correlation function, rather than only the imaginary 
part, is shown to be necessary for accurate formulation of the sum-rules. The charmonium hybrid mass prediction 
is discussed within the context of the XY Z resonances. 

CD. 

^ ■ 1 Introduction 



■ Quantum chromodynamics (QCD) seems to allow for hadrons which contain explicit gluonic degrees of freedom i.e., 
I> ■ glueballs and hybrids. Despite decades of dedicated effort by experimentalists and theoreticians, no such state has 
. been conclusively identified. Within the context of heavy quarkonia, hybrids can make their presence known in two 
] ways: through J^'~^ quantum numbers that are not permissible for conventional quarkonia (so called exotic hybrids), 
and through an overpopulation of states with conventional (non-exotic) quantum numbers. In this work we focus on 
the latter scenario. 

A promising area in which to search for such hybrid states has been provided by the recent population boom in 



■ the charmonium sector above _D£)-threshold [U El El SJ [5l [6] . Since 2002, more than a dozen new resonances have 

■ been discovered, the so-called XYZ resonances, mainly by the Belle and BaBar collaborations; however, few of these 
' particles fit neatly with a conventional charmonium meson interpretation 7 . Not surprisingly, there has been much 

speculation that some of the new states lie outside the constituent quark model. 

In this article, we analyze pseudoscalar (J^'^ — ^) charmonium and bottomonium hybrids using a QCD 
Laplace sum-rules approach. The pioneering calculations for heavy quark hybrids were handled with a constituent 
gluon model [8|. Additional computational approaches (relevant to the pseudoscalar sector) include the flux tube 
model ^ which predicts the lightest charmonium hybrids at 4.1-4.2 GeV as well as lattice QCD [10l|TTl[T2] which 
yields a quenched prediction of 4.01 GeV and unquenched predictions of about 4.2 GeV. To our knowledge, the 
only sum-rules literature concerning heavy quark pseudoscalar hybrids is Refs. [Ml [151 US]. As noted therein, the 
sum-rules that were derived demonstrated instabilities when analyzed. With this paper, we aim to address these 
instabilities and update the sum-rule mass prediction by extending previous work |14[ 115] to include dimension-six 
gluon condensate effects. 

In Section [21 we compute the relevant two-point correlation function. We include leading-order perturbative 
contributions as well as contributions stemming from the dimension-four and dimension-six gluon condensates. Using 
these results, we then derive the needed Laplace sum-rules. In Section [3l we analyze the sum-rules using the single 
narrow resonance model and extract ground state mass predictions. Finally, in Section [4l we comment on our 
charmonium hybrid results and interpret them within the context of current experimental data. 
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2 Laplace Sum-Rules for the Pseudoscalar Heavy Quark Hybrids 

The pseudoscalar ( J^*^ = ^) heavy quark hybrid states ean be studied from the following correlation function [14] 
n^^iq) = I I d^xe^'i-%0\T[j^{x)3M] |0) (1) 



— |Q'^"7'^G'°yQ , — -^^^_lva|^G'^^p , (2) 

where Q denotes a heavy (charm or bottom) quark field. Within ([T|) , the longitudinal part lis is of primary interest 
in this work because it probes the states 

n,.(g) = - .9..) n„(g^) + ^n,(g2) . (3) 

The leading-order perturbative and gluon condensate (aG^) — {ctG'^^G'^^) contributions to the imaginary part 
of lis have previously been calculated [TH [TS], but the resulting sum-rule analysis for the mass was unsta- 
ble [TSl [TB]. We extend these results by calculating the leading-order dimension-six gluon condensate (g^G^) = 
{g^ fabcG^^^G^^^G'^^^) contributions to lis. As will be seen below, the (g^G^) contribution is sufficient to stabilize the 
sum-rule mass prediction. The stabilizing effect of (g^G^) has also been observed for the sum-rule analysis of 
1 heavy quark hybrids [T7] . 

We begin by verifying the leading-order perturbative and (aG^) results [HI [15] for lis. Ref. [18| advocates the 
desirability of an independent confirmation of the Ref. [HI [TS] results; as such we have calculated the full expression 
for lis as opposed to simply reproducing the previously-calculated imaginary part [14[ I15| . 

The leading-order perturbative contribution to lis is represented in Fig. [TJ We use the Tarcer |20j implementation 
of loop-integral recurrence relations and tensor structures [21] [22] to express lis in terms of the small set of basic 
integrals given in Refs. [23l |24l|25]. In £) = 4 -|- 2e dimensions in the MS scheme, the perturbative result is 
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180(z - 1) (4z2 - 21z + 10) 3F2 (1, 1, 1; 3/2, 3; z) 

(4) 

+20z (8z^ + + 29z - lO) 3F2 (1, 1, 2; 5/2, 4; z) 



z 
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where m is the quark mass, and non-physical terms corresponding to dispersion relation subtraction constants have 
been omitted. The quantities a and m are implicitly evaluated at the renormalization scale in the MS scheme. 
Standard conventions for generalized hypergeometric functions have been used (see for example Ref. [26|). The 
analytic structure of lis (i.e., a branch starting at — 4m?) is clearly evident from the hypergeometric functions. 
Analytic continuation of the hypergeometric functions in ^ gives 



Imnr'*(g2) = — ^— V^^Vi (30 - 115z + 166z2 + 8z^ + IGz*) 
1207r-^z^ \ 

\ (5) 

-15(-2-h9z-16z2-|-16z^)log[^/I^+^/I] , z>l. 

Calculating the integral representations for ImllP''''' given in [T31 |TS] we find complete agreement with ([5]) . 

The leading order (aG^) contribution to lis is represented in Fig. [2] We choose to calculate this contribution 
using fixed-point gauge methods (see, e.g., Ref. [27| for examples applying these methods), which have been proven 
to be equivalent to plane- wave techniques for correlation functions of gauge-invariant currents |28j Using the same 
loop-calculation methods as for the perturbative contributions, the (aG^) result is 

n?°(g') = ^-^m'ziAz + 2hF, (1, 1; 5/2; z) , (6) 

where non-physical terms corresponding to dispersion relation subtraction constants have been omitted. The imagi- 
nary part of ([5]) 



Imn^G(g2)^ ^ (l + 2z)^l^, z>l (7) 
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Figure 1: Feynman diagram for the leading-order perturbative contribution to Us- The current 
is represented by the symbol. This and all subsequent Feynman diagrams were created with 
JaxoDraw [19] , 




again agrees with the explicit result of Refs. [Ml [T5] . 

The dimension-six gluon condensate contributions represented by the diagrams in Fig. [3] were not calculated in 
Refs. [TH[T3]. Using the fixed-point gauge and loop calculation techniques described above, we find 



jjGGG/ 2^ ^ (fffG^ 

' ' 3847r2 



3847r2(z - 1)^ - + 16 + - 1)(9 - 21.)] . 



The corresponding imaginary part of ([5]) is 

Imn«««(,^) = ^^^^^^[4z^(z-l)-(4z3-6z^ + 2z-l)], z>l. (9) 

At this stage, we note that only the hypergeometric terms contribute to the imaginary part, but as will be shown 
below, the remaining terms do contribute to the QCD Laplace sum-rules because of the single and double poles at 
z = 1. Thus if only the imaginary parts are calculated (as in Refs. [UJ [TS]), there exists the possibility that the 
resulting sum-rule will be inaccurate. 



Figure 3: Feynman diagram for the leading-order (g'^G^) contribution to lis. Additional diagrams 
related by symmetry are not shown. 

The dispersion relation appropriate to the asymptotic (perturbative) behaviour of Us is 

oo 

n, (q2) ^ n,(o) + Q'n'M + ^Q^Km + ^Q'n-(o) -i- q«1 J dt jt^qT^ ■ (lo) 

to 

'^Implementation of fixed-point gauge methods is trivial because the field strength appears in the current l(2)l. 
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where — ~<f is the Eudidean momentum and p[t) is the hadronic spectral function with physical threshold to- 
Note that the high-energy behaviour of ([5]) ensures convergence of the integral. Direct application of the dispersion 
relation is not possible because lis contains field theoretical divergences that are polynomials in and the associated 
subtraction constants on the right-hand side of (|10p are unknown. A related problem is the contribution of excited 
states and the QCD continuum to the integral of p(t) in pUj) . Enhancement of the lowest- lying resonance contribution 
in hadronic systems requires greater high-energy suppression of this integral. 

The established technique for dealing with these issues is the Laplace sum-rules [29l . A family of Laplace 
sum-rules can be obtained from the dispersion relation (1101) through the Borel transform operator B 



«, q2^„ V{N) \dCf^ 

N/Q^=T 

which has the following useful properties in the construction of the Laplace sum-rules: 

B [ao + aiQ2 + . . . a^Q^™] ^ q , (m finite) (12) 



B 



= T(-l)"t"e-*^ , n = 0, 1, 2,... (n finite). (13) 



The Borel transform is related to the inverse Laplace transform via |31) 



oo 

/ (Q2) = j drF{r)e-Q'^ ^ C [F(r)] \b [f (Q^)] = f{t) ^ £"1 [/ (Q^)] (14) 



/ fiQ")^''^'^dQ" (15) 

h—ioo 

where / (Q^) is analytic to the right of the integration contour in the complex plane. 
The theoretically-determined quantity 

£fe(T) = is[(-l)'=Q2fcn,(Q2)] ^ (16) 

leads to the following family of Laplace sum-rules, after application of B to the dispersion relation (|10l) weighted by 
the appropriate power of Q^: 

00 

£k{T) = ^ J dtt''e-'^p{t) , k>0. (17) 

to 

On the right-hand side of (jl7p . we impose the standard resonance plus continuum model 

p{t) = phad (^t) + e{t-so) imnQCD (t) (18) 

where sq represents the onset of the QCD continuum. The resulting continuum contribution 

£^°"*(s, r, So) = - / exp [-tr] ImnQC°(t) dt (19) 

^ J so 

is then moved to the left-hand side of (IT71) . The total QCD contribution 

£Q^°(r,so)^A(r)-/:rMr,so) (20) 
is then related to the hadronic spectral function 

£^^° (r, so) = - / exp [-tr] p'^-'^it) dt . (21) 

^ J to 

We also note that the tensor decomposition in could have been chosen without the overall factor of l/g^ , as 
is done for axial- vector and vector correlators (for examples of each see Ref. [18j). If this convention had been used, 
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the perturbative calculation would have a l/q^e^ divergence, which must be eliminated by additional weights of Q"^ 
in the Laplace sum-rule (|16p . In other words, knowledge of the divergence structure, which is not revealed in the 
imaginary part, places a bound on the lowest-possible weight k in the Laplace sum- rule (jl6p . 

The Laplace sum-rule can now be calculated using the methods described above (see, e.g., Ref. [22] for detailed 
examples of applying inverse Laplace transform techniques), leading to the following results 



QCD 



(t, So) 



io/4ri 



lim 

?7->0+ 



4m 



so/4m 



imn 



GGG 



Imll^'^ (Am^x)] exp (-4TO^rx) dx 

4m2(g3G3^ 



/:?^°(r,so) 



i+n 
{t, So) ■ 



(4to^x) exp (— 4m^Ta;) dx 



1287r2^ 



exp (— 4m^T) 



(22) 



(23) 



Several clarifying remarks on Eqs. (j221 123p are needed. First, the mass m and strong coupling a are implicitly MS- 
scheme running quantities evaluated at a scale /i. However, one often implements renormalization-group improvement 
by setting fi — 1/\/t [33], which must be done after calculating the r partial derivative in (I23p . Second, the rj limiting 
procedure naturally originates from the inverse Laplace transform approach applied to the full result ([8|), and ensures 
cancellation of integration divergences arising from the l/{z — 1) poles of ((9|). Thus, as mentioned earlier, the last 
term in (|22|) requires knowledge of the full {g^G^) contributions and cannot be obtained solely from ImH^*^*^. 



3 Analysis: Mass Predictions for the Pseudoscalar Heavy Quark Hy- 
brids 

We analyze the QCD Laplace sum-rules using a single narrow narrow resonance model 

iphad(^) ^ j-2g _ ^j2-^ (-24) 

In this approximation, the sum-rules become 

and the ^ hybrid mass M is given by the ratio of the first two Laplace sum-rules 

/.QCD / N 

QCD ■ (26) 

Using fairly general arguments, one can demonstrate that the narrow-width mass estimate would overestimate the 
actual mass when resonance width effects are included |34j . Furthermore, the sq — > oo limit provides an upper bound 
on the ratio (|26p . permitting a very robust upper bound on the ground state mass prediction that is essentially 
independent of the QCD continuum approximation and resonance model. 

Before proceeding with the detailed analysis, the QCD parameters will be specified. It is easy to see that the 
quark mass m sets the basic scale of the mass prediction, so it is the most crucial parameter in our analysis. We 
have chosen to focus on sum-rule estimates of quark masses as they would provide the greatest possibility of a 
self-consistent prediction for the hybrid mass. In particular, the following values encompass the MS quark masses of 
Refs. [35||36ll^|38| 

TOc(At = mc) = mc = (1.28±0.02) GeV, (27) 
mf, (a* = mb) =7nh = (4.17 ±0.02) GeV . (28) 

These values are within the Particle Data Group's recommended ranges [55] . 

Since our calculation is leading order, one-loop MS expressions for the renormalization group evolution of the 
strong coupling and quark masses are appropriate. For the hybrid charmonium analysis, the strong coupling is best 
determined by evolution from the r mass, and in the hybrid bottomonium case by evolution from the Z mass: 

= — rJ^^^^ / . \ ' " = 0-33 ; (29) 

l + H^log(^) 

"(m) = / .X , « {Mz) - 0.118 . (30) 
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23a(Mz) 
12-K 



5 



The numerical values of a {Mr) and a {Mz) are based on the determinations of [10], and we use Particle Data Group 
values of the r and Z masses [39j. The scale dependence of the MS masses can then be expressed to the same 
leading-order as 



/ / N \ 12/25 

rricin) = rUc ._ , 
\a[mc), 

/ / \ \ 12/23 

nibi^i) = nib , . 



(31) 
(32) 



We set [i = l/y/r in our sum-rule analysis [33] , 

For the QCD condensates, we use the following determinations of the QCD condensates from heavy-quark systems 



(g^G^) = (8.2± 1.0)GeV^(aG'2) 
(aG^) = (7.5 ± 2.0) x 10"^ QeV^ . 



(33) 
(34) 

For the (g^G^) contributions to (1^^ we find that 77 = 10^'' is sufficient to evaluate the limit. 

We begin with the analysis of hybrid charmonium. We first establish a window for the Borel parameter r for 
which the sum-rule analysis is considered reliable. Following jSOj, we define the quantities 



fcont {t, So) 



fpow So) 



£Q^"(r,.o)//:?^"(r,.o) 
Cr^ir,oo)/Cr^{r,oo) 

^r'*' (t. So) /-Cq'"'* (r, So) 



(35) 



(36) 



where C^^'^^ includes only the perturbative corrections arising from (j22p . The ratio p5p represents the relative impor- 
tance of the continuum to the sum-rule ratio (|26l) while p6p represents the relative importance of non-perturbative 
(power-law) effects. Consistent with [30 , we define the window of sum-rule validity by fcont > 0.7 (i.e., the continuum 
contribution does not exceed 30%) and 0.9 < fpow < 1-1 (i-e., the non-perturbative contributions do not exceed 10%). 
Figure U show the resulting constraints on the Borel parameter r for the optimum value of sq to be discussed below. 
The resulting region of validity 2.6 GeV^ < 1/t < 4.8 GeV^ is comparable to the window established for the 1^~ 
charmonium hybrid [17]. If we change to the pole scheme for the charm quark (with a pole mass m^"'*^ — 1.71 GeV 
[39j ) ■ the sum-rule window diminishes considerably and thus the sum-rule is less reliable than in the MS quark mass 
scheme, consistent with the findings of [41]. 

fpow ' /cont 



2.5 3.0 3.5 4.0 4.5 5.0 (GeV ) 

Figure 4: The quantities fpow ij, so) (solid line) and fcont (r, so) (dashed line) for hybrid char- 
monium are shown as a function of the Borel scale 1/t for the optimized value so = 23.0 GeV^. 
Central values of the QCD parameters have been used. 

The optimized value of sq and mass prediction is obtained by finding the minimum value sq = 19 GeV^ for which 
the ratio ([26]) stabilizes (in this case, a minimum) at a t value within the so-dependent region of validity. We thereby 
establish a region sq > 19GeV^ and 1.6 GeV < l/y/r < 2.0 GeV for locating an optimized prediction. We then 
search for the optimized mass prediction M and so that minimize the quantitjij 



x'iso)=j:{T7 

j 



/:?""(r„so) 



- 1 



(37) 



A fit based on the quantity 



1 



^QCD 



leads to virtually identical optimizations. 
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where 1.6 GeV < 1/ .jTj < 2.0 GeV. This procedure results in sq 
mass 3.82 GeV. In Fig. [5] we show the optimized ratio for sq - 
including the sq ^ oo limit used for obtaining mass bounds. 



= 23.0 GeV^ and the predicted charmonium hybrid 
23 GeV^ in addition to larger and smaller values, 



VXi/Xo (GeV) 



2.8 3.0 3.2 3.4 3.6 3.8 4.0 (GcV ) 

Figure 5: The ratio (t, so) / C'^'^^ (''", so) for hybrid charmonium is shown as a function of 

the Borel scale 1/r for the optimized value so — 23GeV'^ (solid curve). For comparison the ratio 
is also shown for so — 28GeV'^ (upper dotted curve) and sq = 19GeV^ (lower dotted curve). The 
uppermost dashed curve represents the sq — >■ oo limit corresponding to the bound M < 3.96 GeV. 
Central values of the QCD parameters have been used. 



Uncertainties in the mass prediction resulting from the QCD input parameters are dominated by variations of the 
charm quark mass ^7} and (g^G) while {aC^) variations (|34p are relatively stable. The analysis is also stable 
under an alternative choice of renormalization scale (^ = m'^°^'^ = 1.71 GeV). Adding the uncertainties in quadrature, 
we find the predicted value of the charmonium hybrid mass to be M = (3.82 ± 0.13) GeV. The influence of {g^G) 
on the mass prediction corroborates our key observation that the dimension-six condensate effects are essential for 
stabilizing the mass prediction. The basic scales of our analysis align well with Ref. [T7] which also included effects 
of (g^G^) to find a 1 charmonium hybrid mass of approximately 4.4 GeV for sq « 26 GeV^. By contrast, Ref. P^fJ 
was not able to obtain optimized mass predictions for and 1 hybrid charmonium, so we speculate that the 
dimension-six condensate (g^G^) is a necessary component of sum-rule analyses for heavy quark hybrids. 

For hybrid bottomonium a simple scaling behaviour in moving from the hybrid charmonium to bottomonium 
systems will not occur because in addition to a function of jvr? that would lead to scaling behaviour, there are 
differing pre- factors of the quark mass for each contribution [see Eqs. (|4l6l8p ]. This will lead to differing weights of 
the various perturbative and non-perturbative contributions, and hence there are intrinsic field-theoretical differences 
between hybrid charmonium and bottomonium systems. 

The details of the bottomonium hybrid analysis proceeds in a very similar fashion as the charmonium hybrid case. 
The sum- rule window of validity for the optimized sq = 140 GeV^ is shown in Fig. [51 and using the same methodology 
described above, we establish the region sq > 115 GeV^ and 2.9 GeV < l/y^ < 4.2 GeV for locating an optimized 
prediction. The optimization procedure described above yields so = 140 GeV^ and the predicted bottomonium hybrid 
mass 10.64 GeV shown in Fig. [7] in addition to larger and smaller sq values, including the sq — >■ oo limit used for 
obtaining mass bounds. The bottom quark mass and (g^G^) variations (I28I33|) still dominate uncertainties resulting 
in the final mass prediction M = (10.64 ± 0.19) GeV. Once again, the basic scales of our predictions are in good 
agreement with the 1~~ results [T7] . 




' ' ,5 ' ' ' ' 20 ' ' l/r(GeV^) 

Figure 6: The quantities /pom (t, so) (solid line) and fcont{T,so) (dashed line) for hybrid bot- 
tomonium are shown as a function of the Borel scale 1/r for the optimized value so — 140 GeV^. 
Central values of the QCD parameters have been used. 
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Figure 7: The ratio £j {t,so) / C-q {t, so) for hybrid bottomonium is shown as a function 
of the Borel scale 1/r for the optimized value sq = 140 GeV ^ (solid curve). For comparison the 
ratio is also shown for so = 155 GeV^ (upper dotted curve) and so = 116 GeV^ (lower dotted 
curve). The uppermost dashed curve represents the sq — s- oo limit corresponding to the bound 
M < 10.84 GeV. Central values of the QCD parameters have been used. 



4 Conclusions 

In this paper we have calculated the leading order perturbative, (aG^), and (g^G^) contributions to the pseudoscalar 
(J^*-^ = '") heavy quark hybrid correlation function. A full calculation of the perturbative and (aG^) terms has 
been performed, and the imaginary parts confirm the results of Refs. [131 [H]. However, the (g^G^) contributions 
have not previously been calculated, and the full contribution to the correlation function was needed because the 
imaginary part was not sufficient to determine the Laplace sum-rules. 

In the absence of the {g^G'^) contributions, a stable Laplace sum-rule prediction of the charmonium and bot- 
tomonium hybrids was not achieved even with sophisticated coupled sum-rule methods [15l|T6]. However, the (g^G^) 
effects are able to stabilize the Laplace sum-rule mass analysis and result in the predictions M = (3.82±0.13) GeV for 
the charmonium hybrid and M = (10.64 ±0.19) GeV for hybrid bottomonium. The uncertainties in our mass predic- 
tions only include effects of the QCD input parameters; we make no attempt to estimate the effect of higher-loop or 
other theoretical uncertainties. We emphasize that the (g^G^) uncertainty is clearly observable, demonstrating that 
the dimension-six contributions are significant enough to stabilize the analysis. Ref. [17] previously found a similar 
stabilizing effect of the (g^G^) contributions for mass predictions of 1 hybrid charmonium and bottomonium, and 
the sum-rule scales of the 1 and ^ systems are in qualitative agreement. 

The results of our analysis may have implications concerning the Y(3940) first observed by the Belle Collabora- 
tion [35] and seemingly verified by the BaBar Collaboration [33] although at the significantly lower mass of 3915 MeV. 
There seems to be an emerging consensus that the Y(3940) and the X(3915) [44j are the same particle whereas the 
Y(3940) and the X(3940) gS] are distinct (HHIliH]. In what follows, we adopt this point of view. 

As first noted in the paper announcing its discovery [421 , the Y(3940) is a legitimate charmonium hybrid candidate. 
It is observed in B decays which, as argued in [49) . are thought to be prime charmonium hybrid hunting grounds. 
Also, to date, the only hadronic decay mode detected is to ojj/ijj [371 HH], an observation difficult to reconcile 
with a conventional charmonium meson assignment considering the kinematically- allowed DD and DD channels. 
Such peculiar decay signatures are, however, consistent with a hybrid interpretation as there exists a fiux tube 
model-inspired selection rule which heavily suppresses hybrid decays to pairs of S-wave mesons [SUl jSH (55] . 

The ujJ/ip decay mode allows for a straightforward identification of the Y(3940) as an isosinglet with C = +. 
Unfortunately, the assignment is not so simple; additional effort is required to identify the needed spin and parity 
quantum numbers. 

Optimistically assuming that the Y(3940) will eventually be identified as a pseudoscalar state, comparing our 
mass prediction of 3820 MeV (and 130 MeV uncertainty) to the measured value of 3915 MeV provides additional 
evidence in favour of an interpretation of the Y(3940) as a charmonium hybrid or at least as a resonance admitting 
a significant hybrid component. 
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